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Abstract. The cross-correlation search for gravitational wave, which is known as 'radiometry', 
has been previously applied to map of the gravitational wave stochastic background in the sky 
^ I and also to target on gravitational wave from rotating neutron stars/pulsars. We consider 

. the Virgo cluster where may be appear as 'hot spot' spanning few pixels in the sky in 

radiometry analysis. Our results show that sufficient signal to noise ratio can be accumulated 
with integration times of the order of a year. We also construct numerical simulation of 
radiometry analysis, assuming current constructing/upgrading ground-based detectors. Point 
spread function of the injected sources are confirmed by numerical test. Typical resolution of 
radiometry analysis is a few square degree which corresponds to several thousand pixels of sky 
mapping. 
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^ • 1. Introduction 

Current constructing and upgrading ground-based gravitational wave detectors may reach the 
strain sensitivity in order of 10~^^[1/\/Hz] in near future. These detector sites are dislocated, 
and make possible to have a independent observation of gravitational waves. Since the dislocated 
detectors may not have a correlation in noise component but have a correlation in gravitational 
waveforms, cross-correlation product by two or more detectors make possible to search for 
unknown gravitational waveforms. This is useful property in burst gravitational wave search for 
short duration, or in stochastic gravitational wave search for long duration. Distance between 
detectors also make possible triangulation for signal measurement. Time delay between detectors 
suggest the incident direction of the gravitational wave. It is used in short duration signal 
analysis. On the other hand, an employment of the time delay term in long duration signal 
search or stochastic gravitational wave which come from any directions is not simple. 
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However, we can extract the signal strength if we resolve Earth's rotation and convoluted 
signal from anywhere in the sky. Using multiple detectors and long time integration, the 
deconvolution will be possible. The directional cross-correlation search for gravitational wave is 
proposed as 'radiometry'. 

Radiometry and any other triangulation analysis require long distance among detectors, and 
variation of arrangement (zenith direction, and azimuthal rotation of interferometers) make 
possible to cover the sky and resolve the polarization. Recent starting of the construction of 
LCGT[lJ in addition to previously constructed advanced detector sites of LIG0[2j and Virgo[3] 
encourage the development of radiometry analysis. 

1.1. Radiometry Filter 

Gravitational wave radiometry is a kind of stochastic gravitational wave search, but it is a 
directional analysis. Radiometry analysis assumes that gravitational wave come from particular 
direction as a point like source or a hot spot of sky. The radiometry filter uses at least two 
detector observation data, and do not require the waveform. Gravitational wave coming from 
particular direction will have a time delay of arrival at two detectors shown in FigHJ Since 
gravitational wave signals on both detectors are coherent in same polarization component, we can 
extract gravitational wave signals in the cross-correlation product from two detectors observation 
data with appropriate time delay. Moreover, signal-to-noise ratio for gravitational wave will 
be emphasized with long integration time, because the two detector noises are independent. 
Changing the assumption of source direction step by step, we can search gravitational wave 'sky 
map' by radiometry analysis. 




Figure 1. Ideal setup of a radiometry 



Radiometry analysis is given and studied in previous works 0] [5]. We explain the radiometry 
filter briefly here. The source direction or particular direction of search can be represent as unit 
vector Ct. Two detectors are noted as 1 and 2, and a base-line from detector 1 to detector 2 can 
be given as vector Ax. An arrival time delay for gravitational wave between two detectors is 

, where c is a speed or light. Detector output from each detector S2(t — ) m 

c c 
time domain have same gravitational wave timing. In Fourier domain, this time delay represent 

as a phase difference of two detectors output S2{f) in Fourier transform of si{t) and 52 (t). 

Therefore, we can define the radiometry filter as following. A directed overlap reduction function 

7(/, Q) can be defined to correct phase difference as 
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where F]^2 antenna pattern [6J of each detector 1 or 2 for each two polarization A = + or x 
of gravitational wave. To make optimal filtering, we employ two detectors noise power spectrum 
^1,2 (/) respectively, and assume that gravitational wave spectrum as H(f). Radiometry filter 
can defined as 

and filter output is 

/oo 
df si{frs2if)Q{f,^). (3) 
-oo 

The base-line Ax will change relatively in celestial coordinate according to Earth's rotation 
in the case of ground-based detectors. Therefore, we calculate AS at short time slice ('chunk') 
during t — At/2 and t ^ At/2. Fourier transforms of detector signals can be represented using 
short Fourier transform of si 2{t'J)^ ^^^d Q{f) can be denoted as Q{t^ /). 

/oo 
df ~sl{t-f)h{t-J)Q{tJM). (4) 
-oo 

Processing many chunks for long duration observation data, we can get the statistics of AS. 
Since the gravitational wave component will appear in real part of AS, we employ a mean ^(f^) 

= (5R[A5]) , (5) 

and a standard deviation a{yt) for a real part of AS. Signal-to-noise ratio of particular direction 
can be given as 

p{h) = p{n)/a{Ci). (6) 



2. Hot Spot of Gravitational Wave Radiometry 

Radiometry analysis might find a hot spot in sky map. Hot spot may be formed by concentration 
of unresolved huge number of astronomical sources. Here is one of the promising targets the 
Virgo cluster which could contain ^ 10^^ neutron stars. Our results show that sufficient signal 
to noise can be accumulated with integration times of the order of a year. We published this 
study in the reference [7J. If an average ellipticity of neutron stars are 10~^, Virgo cluster might 
be appear as a hot spot with observation of order of year, using 2nd generation ground-based 
detectors. The numerical estimation of expected GW spectrum and signal-to-noise ratio of Virgo 
cluster hot spot is explained in this reference in detail. 

Here we remark the average detector response for the sky. In eq.([T]), magnitude of two 
detector's combined detector response r{Q,t) is give as 

r{n,t) = F+i(f^,t)F+2(f^,t) + Fxi(f^,t)Fx2(f^,t). (7) 

Since T{(l,t) will change according to Earth's rotation, we have to use average during one 
day : 

-1 /♦Iday 

(r2)iday(fi) = — / r^iCi,t)dt. (8) 

-^Iday Jo 

(r^) will not depend on a right ascension. Fig|2] displays square root of (F^). For Virgo cluster 
whose center is about at declination 12°4'59'^ LCGT and India site (IndIGO[^) is a good 
combination. Variation of combination will cover the whole sky in radiometry analysis. 
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Figure 2. Average of detector response as a function of celestial declination angle. 



3. Numerical Simulation : Source Injection Test 

To confirme the behavior of a hot spot in the radiometry sky map, we construct numerical test 
of radiometry filter. It consist of the following steps. 

1) set short duration of data stream along time axis as ^chunks' : In real observa- 
tional data processing, this step and following steps 2) 3) 4) are corresponding to split 
time series data into short duration as 'chunks'. Following steps 2)-5) are processed for 
each chunk. 

2) generate gravitational wave detector noise : Assuming design noise spectrum [9] 
[To] [IT] [T2] of each detector and injecting random gaussian noise, we generate time series 
of noise n{t). In this paper, time series data's sampling frequency is 4000 Hz, and chunk 
length is 65.536 sec. In the study of point spread function in next subsection, we skip noise 
generation to use gravitational wave signal only. 

3) generate gravitational wave signal : Gravitational wave signal at source is generated. 
We tested with sinusoidal wave, and random phase waveform which has power spectrum in 
certain frequency band, e.g. flat power spectrum in 300ib30 Hz. 

4) inject gravitational wave signal into detector noises in time domain : 

Gravitational wave signal will incident detectors with Doppler effect, antenna patterns and 
time delay between two detectors. Doppler effect is employed as a frequency modulus of 
gravitational wave on the detector according to the orbital motion and Earth's rotation. 
For incident direction on detectors, antenna patterns F+,Fx are calculated for the signal. 
Also the arrival time difference of the signal on two (or more) detectors are given in this 
step. These condition will change as time goes by according to Earth's rotation. 

5) calculate radiometry filter output AS{t,Q) : We calculate radiometry filter output of 
this chunk in eq.Q using two detectors for the particular direction of a ^pixeV. Changing 

step by step, sweeping all sky, we calculate AS'(t,f^) for all pixels. 

6) repeat step 1) — 5) for next chunk : We iterate steps above for all chunks. 



Finnaly, 

7) sum up AS and calculate signal-to-noise ratio p{(V) of each pixel : We sum up 

each cunk's output AS'(t,fi) with respect to each pixel. Then we get ji{Ct)^(j{(V) and 



calculate signal-to-noise ratio p{(l) in eq.®. p{Ct) are displayed as pixel of sky map with 
colorized according to its magnitude as Fig|3l 



FigH] displays a schematic data flow of our numerical test of radiometry filter of above steps. 
We tried possible combinations with current detectors, but FigH] displays LCGT and advanced 
LIGO for example. For studying point spread behavior in following sections of this paper, we 
display also in the case of LCGT and advanced LIGO at Livingston. 
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Figure 3. Example of sky map of radiometry 
output. This map is displayed in celestial 
coordinate. 



Figure 4. Data flow of numerical test of 
radiometry filter 



3.1. Point Spread Function 

How spread the point source on the sky is important to evaluate the feasibility of searches. It is 
called as ^point spread function'. Point spread function depends on gravitational wave frequency 
or spectrum, a position on the sky (right ascension, celestial declination), and for detector 
orientation on the Earth. In point spread function study in this subsection, we used noiseless 
detector output, but other filter setting is same to explanation above. Thus, the magnitude 
of the gravitational wave in this section is arbitrary scale. Test injection of gravitational wave 
signal has a flat spectrum around /o [Hz] as 



for average spectrum. Injected signal for each chunk fluctuates as stochastic signal whose phase 
is random. 

Fig|5] shows how point spreading for some frequency samples in single chunk process. In single 
chunk, we can find of wide and multiple spreading image from the injection of point source, since 
the different time lag makes aliases. Spatial interval is wide for low frequency band, and narrow 
for higher frequency band. But in general, it is hard to identify the exact position of the point 
source in single chunk map, even if the gravitational wave is so strong. 

With Earth's rotation, aliases will change its orientation since detector position will change 
relatively to the source. Therefore, during long integration time, aliases will be reduced and true 
position will grow up with stacking filter output coherently. FigE] displays radiometry map with 
one day integration for closed up narrow region around injection sources. The figure displays 
variation in case of declination angle of 5 = 15°,45°and 75°. The point spread function is not 
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Figure 5. Radiometry filter output of single chunk with various frequency signal. Star symbol 
shows true injection position in the sky. Gravitational wave has a flat spectrum around fo [Hz] 
with zb 10% of fo [Hz] width; e.g. /o=300Hz means from 270Hz to 330Hz. 



depend on right ascension with more than one day integration. According to the combined 
detector response as eq.([2|) with LCGT and LIGO detectors, mid dechnation angle have a 
circular image of spreading. On the other hand, in low and high declination area, point spreads 
as ellipsoid. 




Figure 6. Point spread image in radiometry map (closed up narrow region) with one day 
integration with variation of declination angle. Frequency fo of injected signal is 300 Hz. 
Declination angles of injection point is 5 = 15°,45°and 75° respectively. 



Spread size of the point source is strongly depend on the gravitational wave frequency. Its 

order is determined by a diffraction limit y , where wave length A of gravitational wave and 

\Ax\ 

distance between two detectors |Ax|. Fig|7] displays some examples of spread image of injection 



point at 5 = 60°. The signal of 30 Hz spreads widely. IkHz signal has aliases as concentric 
circle, but the true injected point appear as intense peak. 
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Figure 7. Point spread image in radiometry map with one day integration with variation of 
gravitational wave frequency. Injection point is at 5 = 60°. Central frequency of each figure is 
30Hz, lOOHz, 300Hz and IkHz respectively. 



To evaluate these spreading, we define the spread size using solid angle Afi [sr] as a solid 
angle of area which lager than half hight of the local peak. FigjH] displays relation between Afi 

and the frequency of gravitational waves. Solid line is a y and dashed line is a empirical 

|Ax| 

approximation for higher 6 region : 

f 1 c Y f80-S [deg.]\ , , 
^ =^ X ^7-^ sr . (10) 

In the FigiHl we also display corresponding apparent diameter for ^-axis in case of circular 
image. 

3.2. Optimal Pixel Size 

Point spreading suggests on optimal size of pixels. It is possible to employ fine pixels as similar 
size of point spread function, because this is an angular resolution of radiometry analysis. 
However, there is no need to divide the celestial sphere so fine as one pixel size is far smaller 
than point spread function. Figl9] displays optimal number of pixels in the sky when we employ 
as point spread size (solid angle) as pixel size. In this study with LCGT and LIGO (Livingston) 
base line can employ up to several thousand pixels for detectors' good sensitivity frequency band 
around 100 - 300 Hz. 




Figure 8. Point spread size of injected gravitational waves. Points are injected gravitational 
wave signal with different central frequency/o and variation of declination angle 6. x-axis is a 
gravitational wave frequency, and ?/-axis is a solid angle of point spreading, ^-axis is also displays 
corresponding apparent diameter in case of circular image. 




Figure 9. Optimal number of pixels in the sky when we employ as point spread size (solid 
angle) as pixel size. Solid line is corresponding to the diffraction limit size. 



4. Summary 

We studied radiometry analysis for stochastic gravitational waves from astronomical origin, 
which might have an anisotropy or local excess. Virgo cluster has a possibility to appear as a 
'hot spot' on the radiometry map of the sky with advanced and third generation gravitational 
wave detectors. We also studied numerical simulation of radiometry filter, and displayed the 
spreading image for point like gravitational wave sources, with LCGT and advanced LIGO 
combination. Also other combinations in global network detectors will be useful to emphasize 
the signal behavior. We are planning to proceed the injection test of our numerical simulation 
studies for multiple sources and area spreading sources. 
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